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Inclusive electron scattering off the deuteron has been measured to extract the deuteron struc- 
ture function F2 with the CEBAF Large Acceptance Spectrometer (CLAS) at the Thomas Jeffer- 
son National Accelerator Facility. The measurement covers the entire resonance region from the 
quasi-elastic peak up to the invariant mass of the final-state hadronic system W ~ 2.7 GeV with 
four- momentum transfers from 0.4 to 6 (GeV/c)^. These data are complementary to previous 
measurements of the proton structure function F2 and cover a similar two-dimensional region of 
and Bjorken variable x. Determination of the deuteron F2 over a large x interval including the 
quasi-elastic peak as a function of Q^, together with the other world data, permit a direct eval- 
uation of the structure function moments for the first time. By fitting the evolution of these 
moments with an OPE-based twist expansion we have obtained a separation of the leading twist 
and higher twist terms. The observed behaviour of the higher twist contribution suggests a 
partial cancellation of different higher twists entering into the expansion with opposite signs. This 
cancellation, found also in the proton moments, is a manifestation of the "duality" phenomenon in 
the F2 structure function. 
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I. INTRODUCTION 

Inclusive lepton scattering off the deuteron has pro- 
vided a wealth of information about internal nucleon 
structure and nuclear phenomena. Since a free neutron 
target does not exist, the deuteron is the simplest tar- 
get for the study of neutron structure functions. The 
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weak coupling between the two nucleons in the deuteron, 
corresponding to a large space-time separation, suggests 
that the nucleus can be described as a non-relativistic 
proton and neutron moving inside some mean potential. 
To this end, the momentum distribution of the deuteron 
was established with high precision from the quasi-elastic 
and NN reactions, and approaches describing the Fermi 
motion of the nucleons were developed. This naive pic- 
ture, however, was superseded by experiments when the 
European Muon Collaboration (EMC) discovered devia- 
tions of the measured nuclear structure function F2 from 
that of a free proton and neutron convoluted with Fermi 
smearing (a phenomenon known as the EMC-effect Q]). 
Different attempts to explain the EMC-effect have been 
undertaken, but without reaching a unified and, there- 
fore, definitive description. 

A Quantum Chromodynamics (QCD)-based approach, 
on the other hand, can handle the nuclear structure func- 
tions in a model-independent manner. This method is 
based on the Operator Product Expansion (OPE) of 
the structure function moments, whose Q^-evolution is 
known in QCD at leading twist and some fixed order in 
as- Note that QCD predictions on the Q^-evolution are 



3 



target-independent and nuclear effects are only to modify 
the normalization of moments, known in few cases from 
Lattice QCD simulations. However, the leading twist 
picture is valid only in Deep Inelastic Scattering (DIS) 
and for not very large x values. Beyond these kinematic 
boundaries, e.g. in the region of large x and moderate Q^, 
new poorly established physics appears. This kinematic 
domain has a particular interest because multi-parton 
correlation phenomena manifest themselves as deviations 
from perturbative QCD (pQCD) predictions. For a nu- 
clear target some part of these deviations will be related 
to the nucleon interaction with its environment. 

In the one photon exchange approximation, the cross 
section for inclusive electron scattering off a nucleus is 
described by the total absorption of the virtual photon 
by the nucleus. The optical theorem relates the total vir- 
tual photoabsorption cross section to the forward Comp- 
ton scattering amplitude of the virtual photon on the 
nucleus. The latter amplitude, in general, can be repre- 
sented as a product of two hadronic currents separated 
by a certain space-time interval C^- In the Bjorken limit, 
the interval ^ (while the light-cone component is 
fixed) allows one to apply the OPE to the product of non- 
local hadronic currents. This leads to a relation where 
the moments M^^ (Cornwall-Norton definition 0) of 
the nucleon (nucleus with mass A) structure functions, 
defined as: 



M^^(g2)=/ dxa;("-2)i^2(x,Q2), „ > 2, 



n even, 
(1) 

are expanded as a series of inverse powers (twists) of the 
four-momentum transfer (for details see Ref. p]). A 
study of the dependence of these moments, therefore, 
would permit one to isolate different terms of this series, 
each representing distinct physical processes in QCD. 
The first term in the expansion represents the leading 
twist, i.e. the limit of asymptotic freedom, while higher 
terms account for the interactions among partons inside 
the nucleon. The contribution of these multi-parton cor- 
relations to the nucleon wave function increases in the 
region of large x (corresponding to high moment order 
n) and low Q^. 

For the proton, a careful study of the multi-parton 
correlation contribution, which included a global analy- 
sis of all the world's data on the proton structure function 
F2, has been recently described in Ref. Q. This analysis 
could not be done for the deuteron structure function be- 
cause of the lack of data in the region of the quasi-elastic 
peak, and because of the scarcity of data in the reso- 
nance region. However, a previous analysis from Ref. 0, 
based on fits of the structure function F2, showed a mod- 
ification of the scaling behaviour of the nucleon structure 
function F2 in the nuclear medium. The Hall C Collabo- 
ration at Jefferson Lab has recently provided high quality 
data in this kinematic region 0| , but the exclusion of the 
quasi-elastic peak in this measurement prevented further 
studies of these data in terms of QCD. 

In this paper we report on a measurement of unpolar- 
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FIG. 1: Experimental data on the deuteron structure 
function F2(x,Q^) used for the moment evaluation in the 
CLAS kinematic region. The points show world data from 

Refs. II IIS a 113 111 El El El El El El El El 113 123. 

The shaded area shows CLAS data. 



ized inclusive electron scattering from deuterium taken 
with the CLAS detector in Hall B at Jefferson Lab. The 
data span a wide continuous two-dimensional region in 
X and (see Fig. EJ. The F2 structure function of the 
deuteron was extracted over the entire resonance region 
{W <2.7 GeV) below Q'^ = & (GeV/c)^. This measure- 
ment, together with existing world data, allowed for the 
first time the evaluation of the first four F2 moments of 
the deuteron down to 0.4 (GeV/c)^. 

In section HTl we review the F2 moments in the frame- 
work of pQCD. In section lTTll we discuss improvements in 
the data analysis implemented since the first unpolarized 
inclusive measurement at CLAS, along with some details 
of the evaluation of the moments. For other details of 
the analysis we refer to Ref. 0. Finally, in Section IIVI 
we discuss the interpretation of the results. 



II. MOMENTS OF THE STRUCTURE 
FUNCTION F2 

Measured structure functions for a free nucleon target 
in the DIS regime are related to parton momentum distri- 
butions of the nucleon. For a nuclear target this relation 
is not direct since it is necessary to account for effects of 
the Fermi motion, meson exchange currents, off-shellness 
of the nucleon and final state interactions (FSI) . Never- 
theless, the OPE of the structure function moments of 
the nuclear structure function is still applicable in the 
same way as for the free nucleon. The n-th Cornwall- 
Norton moment of the (asymptotic) structure func- 
tion F2{x, Q^) for a massless target can be expanded as: 
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_ (2) 

where fc = l,2,...,oo, /i/ (/^r) is the factorization (renor- 
malization) scale^, ^ is an arbitrary reference scale, 
OnrifJ'r) IS the reduccd matrix element of the local oper- 
ators with definite spin n and twist r (dimension minus 
spin) which is related to the non-perturbative structure 
of the target. EnrilJ^n l^f ,Q'^) is a dimensionless coeffi- 
cient function describing the small distance behaviour, 
which can be perturbatively expressed as a power ex- 
pansion of the running coupling constant as(Cp). More- 
over, the leading twist (t = 2) dependence remains 
unchanged with respect to the free nucleon target and 
all the nuclear effects appear either in the higher twist 
terms (r > 2) or in the reduced matrix element OniifJ-r), 
which does not depend on Q^. The non-zero mass of 
the target^ leads to additional AP /Q^ power correc- 
tions (kinematic twists) which mix operators of different 
spin. These target mass corrections can be accounted for 
by use of Nachtmann [2^ moments (Q^) instead of 
the usual (massless) Cornwall-Norton moments. In the 
Bjorken limit /Q"^ terms become negligible and both 
definitions coincide. The Nachtmann moments for the 
deuteron structure function are defined as follows: 

Ja 

3 + 3{n + l)r + n{n + 2)r^ 

(n + 2)(n + 3) ' ^ ' 



where r = -y/l + 4A4^x^ /Q'^ , M is the proton mass and 
e = 2a;/(I + r). 

The evolution of the leading twist term is known for the 
first four moments up to Next-to-Next-to-Leading Order 
(NNLO). However, if one wants to extend the analysis 
down to « and to large where the rest of the 
perturbative series becomes significant, one needs to ac- 
count for additional log arithmic corrections due to soft 
gluon radiation [23, These corrections re-summed 

in the moment space to all orders of as appear due to 
an imbalance of the virtual and real gluon emission at 
a; — > 1. Since the Q^-evolution of the higher twist terms, 
related to quark-quark and quark-gluon correlations, is 
unknown, their logarithmic QCD behaviour is parame- 
terized and the corresponding anomalous dimensions are 
extracted from the data. 

Measurement of the Nachtmann moments [Q"^) in 
the intermediate Q"^ range (0.5 < < 10 (GeV/c)^) 



[1] We are working in the Soft Gluon Rc-summation (SGR) 

scheme |23| . where /i^ = /x^ = . 
[2] In the leading twist approximation the target is a nucleon inside 

the deuteron. 



allows a model-independent separation of the total higher 
twist contribution from the leading twist. Comparison 
of the higher twist contribution in the deuteron to that 
in the free proton provides important insight into the 
nucleon structure modifications inside nuclear matter. 



III. DATA ANALYSIS 

The data were collected at Jefferson Lab in Hall B with 
the CLAS using a liquid-deuterium target with thick- 
ness 0.81 g/cm^ during the electron beam running peri- 
ods in March- April 2000 and January-March 2002. The 
average beam-target luminosity for these periods was 
6x10"^'^ cm~^s~^. To maximize the interval in and x, 
data were taken at two different electron beam energies: 
= 2.474 and 5.770 GeV. The accumulated statistics 
at the two energies is large enough (> 10^ triggers) to al- 
low for the extraction of the inclusive cross section with 
a small statistical uncertainty (< 5%) in small x and 
bins (Ax =0.009, AQ2 =o.05 (GeV/c)^). 

The CLAS is a magnetic spectrometer based on 
a six-coil torus magnet whose field is primarily oriented 
along the azimuthal direction. The sectors, located be- 
tween the magnet coils, are individually instrumented 
to form six independent magnetic spectrometers. The 
particle detection system includes drift chambers (DC) 
for track reconstruction 26], scintillation counters (SC) 
for time of flight measurements 1271 , Cherenkov counters 
(CC) for electron identification [23] , and electromagnetic 
calorimeters (EC) to measure neutrals and to improve 
electron-pion separation ^2^. The EC detectors, which 
have a granularity defined by triangular cells in a plane 
perpendicular to the incoming particles, are used to study 
the shape of the electromagnetic shower and are longitu- 
dinally divided into two parts with the inner part acting 
as a pre-shower. Charged particles can be detected and 
identified for momenta down to 0.2 (GeV/c) and for polar 
angles between 8° and 142°. The CLAS superconduct- 
ing coils limit the acceptance for charged hadrons from 
about 80% a.i 9 = 90° to about 50% at forward angles 
{9 = 20°). The total angular acceptance for electrons is 
about 1.5 sr. Electron momentum resolution is a func- 
tion of the scattered electron angle and varies from 0.5% 
for e < 30° up to 1-2% for > 30°. The angular resolu- 
tion is approximately constant, approaching 1 mrad for 
polar and 4 mrad for azimuthal angles: the resolution for 
the momentum transfer ranges therefore from 0.2 up to 
0.5 %. The scattered electron missing mass [W) reso- 
lution was estimated to be 2.5 MeV for a beam energy 
less than 3 GeV and about 7 MeV for larger energies. To 
study all possible multi-particle states, the acquisition 
trigger was configured to require at least one electron 
candidate in any of the sectors, where an electron candi- 
date was defined as the coincidence of a signal in the EC 
and Cherenkov modules for any one of the sectors. 

The data analysis procedure has been described in de- 
tail in Refs. 0,113. Therefore, in this article we focus 
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on changes and improvements in the analysis. The most 
important improvements, leading to a significant reduc- 
tion of the estimated systematic uncertainties relative to 
those of Ref. 0, are described in the following sections. 



A. Electron Identification 

The pion contamination observed in Ref. ^ in the elec- 
tron candidate sample was found to be due to random 
coincidences between a pion track measured in the DC 
and a noise pulse in a CC photomultiplier tube (PMT) 
(typically corresponding to one photo-electron). These 
coincidences can be greatly reduced by means of better 
matching between the CC hits and the measured tracks. 

Each CLAS sector consists of 18 CC segments, each 
containing two PMTs. Therefore, the probability of a 
coincidence is the product of probabilities to have a noise 
signal in one of 36 PMTs together with a negative pion 
track in a time interval, At =150 ns, which corresponds 
to the trigger window time. The average CC PMT noise 
rate R^^^^ in CLAS was measured to be «42 kHz. For 
our typical running conditions, the average rate of nega- 
tively charged hadrons within our geometrical EC fiducial 
cuts (used to ensure the shower is fully contained with 
the detector) that have an appropriate EC signal is of 
the order of w2.3 kHz. This gives an estimate of 
the possible contamination: 

j^lphe ^ JlPMTj^h-^^ w 15 Hz , (4) 

which should be compared to the electron rate Ri250 
Hz, using the same cuts. Therefore, the expected con- 
tamination is of the order of 6% overall. In contrast, for 
small momenta p < 1 (GeV/c) and large scattering an- 
gles, e > 30°, R^- K.1.7 kHz and R^- «100 Hz, resulting 
in a contamination of 12%. 

In order to reduce the contamination of the coinci- 
dences between a hadron track and a CC noise signal, 
we applied geometrical and time matching requirements 
between the CC signal and the measured track in the 
following way: 

1. we defined a CC projective plane, an imaginary 
plane behind the CC detector where Cherenkov 
radiation would arrive if it were to propagate the 
same distance from the emission point to the PMT 
without any reflections in the mirror system; 

2. for each CC segment we found the polar angle from 
the CLAS center to the image of the CC segment 
center and to the images of the CC edges; 

3. the impact point and the direction of the track in 
the SC plane, as measured in the DC, were used to 
obtain the measured polar angle 9 in the projective 
plane for each electron candidate event. We fitted 
the 9 distributions separately for each CC segment 
in order to extract their measured width (jp; 



4. for each CC segment we applied a cut: 

\Otrack — 9hit \ < 3crp (5) 

which was intended to remove those electron event 
candidates for which the track impact point in the 
CC was far away from the segment where the hit 
was detected. In Fig. |21 an example of the 9 dis- 
tribution is shown for one segment and the cut ap- 
plied is indicated by dashed lines. To clearly iden- 
tify the contribution of the coincidences and check 
the efficiency of the cut we separated events out- 
side the single photo-electron peak (which contains 
most of the pion contamination) by applying a cut, 
Nphe > 2.5. These electron candidate events with 
reduced pion contamination are shown in Fig.|21by 
the hatched histogram. The difference between the 
empty and the hatched histogram in the figure is 
therefore mostly due to the pion contamination; 

5. to perform time matching between the CC and SC 
hits of the electron candidate event we studied the 
distribution of the time-offset of the CC signal with 
respect to the SC. For each CC segment we mea- 
sured At^'^~'^'^ , defined as follows: 

iSC iCC 

^tSC-CC ^ ^SC _tCC^ l J ^ (g) 

where t^"-^ and t'~^'~^ are hit times, /'^'^ and l^^ are 
the path lengths from the CLAS center to the hit 
points in the SC and CC, and (3 is the track velocity 
measured in the SC. The distribution of SC-CC 
time-offsets is shown in Fig. The electron-rich 
events outside the single photo-electron peak are 
emphasized by the hatched histogram, and the cut 
applied is shown by the dashed line. The presence 
of the double peak near At^'^~'-^'^ is expected 
due to a small time offset between the two PMTs 
within the CC segment; 

6. The measured photo-electron distribution for the 
electron candidate events was compared to the one 
obtained after all the cuts described above had been 
applied. This comparison is shown by different his- 
tograms in Fig. ^ 

After applying the matching procedure, the pion contam- 
ination was reduced from 30% to 5% in the worst case. 
The remaining contamination is due to events where the 
pion impact point in the CC is very close to the PMT 
with the noise signal. This contribution was removed by 
the same procedure as in Ref. 3 • In order to estimate the 
systematic uncertainties of this correction we compared 
the method to another approach, namely, requiring more 
than 2.5 photo-electrons in the CC and estimating the 
number of missing electrons by an extrapolation, using 
the empirical function from Ref. {3ll | adjusted for each 
particular run period. The difference between the two 
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methods provides an estimate of the systematic uncer- 
tainties in the pion rejection corrections. The total rela- 
tive systematic uncertainty of this correction is kinematic 
dependent; at larger the contribution of pions and the 
corresponding systematic uncertainty are larger. 

To improve electron identification at large momenta, 
the energy deposited in the EC was used. Starting from 
momenta of 2.7 (GeV/c), pions begin to emit Cherenkov 
light in the CC. In this region, however, the EC becomes 
very efficient in the separation of electrons from pions 
(see Ref. ^]), and moreover, the ratio of pion to elec- 
tron yields drops down very quickly. Electrons passing 
through the EC release 30% of their energy on average, 
while pion losses are constant. This EC property was ex- 
ploited at large particle momenta for electron-pion sep- 
aration by selecting particles with the energy fraction 
released in the EC above 20%. More details on this pro- 
cedure can be found in Ref. 4]. Furthermore, pions just 
above the Cherenkov threshold emit a small amount of 
light with respect to electrons of the same momenta, and 
therefore, were also removed by the cut on the number 
of photo-electrons. 
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FIG. 3: The time difference between hits in the SC and the 
CC assigned to the electron candidate track, corrected for the 
distance traveled from the SC to the CC. The meaning of the 
two histograms is the same as in Fig.|21 




FIG. 2: The difference in the projective polar angle Q between 
the hit position in the CC and the impact point of the electron 
candidate track. The hatched area shows events with reduced 
pion contamination in the electron candidate sample. The 
reduction of the pion contamination is obtained through an 
additional cut: the number of photo-electrons in the CC > 
2.5. 



B. e e pair production 
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FIG. 4: The spectrum of photo-electrons measured in the CC. 
The hatched area represents the CC spectrum after applying 
the matching cuts described in the text. 



surements were taken with an out-bending torus field'^. 
This provides the possibility to extract the contribution 
of e~ pair production directly from the data. To do 
so, we demanded the first particle in each event to be a 
positron (positron trigger), i.e. to have positive charge 
and hits in both the CC and the EC. We applied ex- 
actly the same cuts and corrections described above and 



The most important source of e^e^ pairs in the CLAS 
is TT*^ production followed cither by Dalitz decay to 7e^e^ 
or by 77 decay with one of the photons converting to 

e~ . For the data set at higher beam energy some mea- 



[3] Normal setting of the CLAS torus magnet bends electrons in the 
forward direction along the beam (i.e. in-bending). The inverse 
magnetic field configuration is called out-bending. 
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in Ref. |j| to the positron trigger data. Following the 
procedure described in Ref. |32| an additional severe cut 
on the number of photo-electrons measured in the CC 
was applied [Nphe > 4) to both the electron and the 
positron trigger rates. The ratio /e~ obtained in this 
way is shown in Fig. in comparison with calculations 
described below. For the higher beam energy data set 
we therefore subtracted the measured e'^e~ background 
accounting for the statistical uncertainties only. 

Since positron data are not available for the lower 
energy data set, the pair production background pro- 
cesses were estimated according to a model developed 
by P. Bosted and described in Ref. [s^. Bosted devel- 
oped a computer code based on the Wiser fit of inclu- 
sive pion production. This model was carefully checked 
against measurements of the inclusive e^/e~ ratio in dif- 
ferent CLAS runs on polarized proton and deuteron tar- 
gets 32] , and it appeared to be in good agreement (within 
30% relative uncertainty) with the measured ratio. The 
value of the correction is assumed to be equal to the ra- 
tio of the inclusive production cross section over the 
fit of the deuteron inclusive cross section c^f^, including 
radiative processes (the tail from the elastic and quasi- 
elastic peaks, bremsstrahlung and the Schwinger correc- 
tion). This correction factor is given by: 



(7) 

where (Te+ is the inclusive production cross section 
and (T^^ is the fit folded with radiative processes. The 
index "M" here refers to the model cross section used in 
the event generator. 

To estimate systematic uncertainties in the calcula- 
tions, we compared the calculated e+e~-pair production 
contribution to the measured one for the higher beam 
energy data. The difference, shown in Fig. [S] has been 
parameterized as a function oiy = v/E and is given by: 



(8) 



where E is the beam energy and v is the energy of the 
virtual photon in the Lab frame. The systematic un- 
certainty estimated for the lower beam energy does not 
exceed 3.5%. 



C. Momentum corrections 

A new correction procedure has been developed in 
Ref. [Ill . This procedure is based on studying completely 
exclusive reactions with overdetermined kinematics in or- 
der to correct simultaneously both momenta and polar 
angles 9 for all particles (azimuthal angle 4> was assumed 
to be correct). The momentum change AP were param- 




FIG. 5; The contribution of e^e pair production events 
in the inclusive cross section at = 1.775 (GeV/c)^ The 
points show the measured quantity 1 — e^/e~, which repre- 
sents the number of electrons inelastically scattered off deu- 
terium to the total number of measured electrons. The curve 
represents the calculations from Ref. |33[ . 



eterized as: 

^^\{E + F<i))^^^ + {G + H(t))sine\ (9) 
F cos © 



P 



{Jcose + Ksm9 + L sin 26) 



qBrorus 

iMco80 + Nsm0 + Osm20)4> , 

where q is the particle charge (in electron charge units) 
and 



Btotus = 0.76- 



Torus 



sin^ 46* 



33756* 



< 7r/8 (10) 



Torus 



0.76 



Ito 



337561 



> vr/ 



Itotus is the Torus magnet current in A and the polar 
angle change is given by 



Ae* = (A + Bq 



COSf 



+ {C + D(t>)smB (11) 



14 free parameters {A — O) were determined separately 
for each CLAS sector from fitting the sum of the squared 
differences of the initial and final four-momentum com- 
ponents in the reactions: 



e -f p ^ e +p 
e+p-^e +p + TT TT 
e + d ^ e' + p' + p" + IT' 



(12) 



During the fitting, the hadron momenta were corrected 
for the energy losses inside the tar get, and a better beam 
energy determination was applied |3J|. 
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Systematic errors of the momentum correction proce- 
dure were studied by means of a comparison between 
data and simulations. Simulated spectra were shifted as- 
suming a constant AP shift that equalized the positions 
of the quasi-elastic peaks. The maximum deviation in 
both sets is of the order of 5 MeV, which is smaller than 
the CLAS resolution. Also the position of the neutron 
peak in the e'n^ missing mass spectrum was checked and 
found to be in agreement with the nominal neutron mass 
value within 3 MeV. The resulting difference in terms 
of the F2 structure function between the direct simula- 
tion spectra and the shifted ones was assumed to be an 
estimate of the corresponding systematic error. 



D. Simulations 

Determination of the acceptance and efficiency correc- 
tions was based entirely on the Monte Carlo (MC) sim- 
ulations developed for CLAS. Moreover, the systematic 
uncertainties of these corrections were estimated from a 
comparison of MC simulations with experimental data 
using a realistic model in the event generator. In short, 
the procedure is the following: we generated events with 
the event generator describing elastic, quasi-elastic and 
inelastic eZ3-scattering processes including radiative cor- 
rections; these events then were processed with GEANT 
based CLAS software simulating the detector response; 
after that the standard CLAS event reconstruction proce- 
dure was applied to obtained detector signals; finally the 
ratio of reconstructed events to the number of generated 
ones gave a combined efficiency/acceptance correction in 
each kinematic bin. 

The simulations of the detector response were per- 
formed in the same way as described in Ref . . The fol- 
lowing improvements and changes for electron-deuteron 
scattering were implemented: 

1. Electron scattering events were generated by a 
random event generator with the probability dis- 
tributed according to cr^^, described in Ap- 
pendix. The values for the elastic and inelas- 
tic cross sections for electron-deuteron scattering 
were taken from existing fits of world data, in ref- 
erences "SS^ and respectively. The contribution 
from internal radiative processes was added accord- 
ing to calculations [s^. 

2. The event rate obtained in the simulations was 
then compared to the data, preserving the original 
normalizations (accumulated Faraday Cup (FC) 
charge for the data and the number of generated 
events over the integrated cross section of the event 
generator for simulations). These normalized yields 
do not include acceptance, efficiency and radia- 
tive corrections. The simulated events passed all 
cuts: the fiducial cuts, calorimeter cut, event status 
cut (see Ref. for details) and CC matching cut. 
However, e~^e~ pair production and empty target 



backgrounds were subtracted from the data. The 
normalized yield obtained with the same set of cuts 
from the data and simulations were compared and 
found to be in good agreement within w 10 — 15% 
(see Fig. 01 , which is at the level of reliability for 
our cross section models. As one can see below in 
this section, a « 10 — 15% variation of the cross sec- 
tion model in the event generator yields only « 1% 
uncertainty in the final cross section. 

In order to check the absolute normalization of the 
inclusive e^D events, we used elastic e~ D scattering 
data. In the range of considered in this analysis 
(0.47-0.63 (GeV/c)^) , the e~D elastic cross sections 
are known within 5% (see Ref. [stI ) . Recent measure- 
ments performed in Jefferson Lab jH, provide the 
most precise data in this kinematic region. We per- 
formed the simulations of this reaction using the pa- 
rameterization from Ref. [35l | and taking into account 
radiative corrections. The normalized event yield ob- 
tained from simulations (^^j,^) was compared to the 
measured one (45 ) at the lowest values. The 

^ all exp' ^ 

yields, shown in Fig. [71 integrated over the peak, are 
in good agreement within statistical and systematic un- 
certainties. The distortion of the measured peak is due 
to the Wd = ^/M"^ -f 2Mdv — Q^-dependence of the ra- 
diative corrections, which were not taken into account in 
the simulations. This W^i-dependence is generated by 
the radiative tail from the e~ D elastic peak which in our 
approach to radiative corrections begins at higher Wd 
due to the cut on soft photons. The efficiency obtained 
from simulations and checked against the elastic scat- 
tering data appear to be approximately constant (about 
90-95%) inside the fiducial region of the detector defined 
by the fiducial cuts. 

The elastic scattering normalized yield was evaluated 
as the number of e~ D coincidences measured in CLAS 
in the Wij-interval from 1.75 to 2 GeV multiplied by the 
corresponding luminosity: 



Lii £ exv.sim 



tot J 1.75 



(13) 

represent the corresponding numbers of 
events (for the measured cross section the empty target 
events were subtracted). In order to clean up the elastic 
data sample, we applied additional cuts on the deuteron 



where N, 



exp,st7n 



identification \M 



2 {exp) 
D 



Ml\< 0.5 GeV2 and the kine- 



matic correlations between the electron and deuteron: 



\\<l)e-4 

COS 9d 



d\ ~ 180° I < 3° and 

{Eo - Pe cos Be 



(14) 



Pd 



< 0.1 , 



where Af^^*^^^^ is the mass of the deuteron measured in 
the SC; Mo is the nominal deuteron mass; 4>e (Oe) and 0£i 
{0d) represent the azimuthal (polar) angles (in degrees) 




FIG. 6; The normalized event yield obtained from the data 
(filled triangles) and simulations (open circles) ai E — 5.77 
GeV, = 2.425 {GeV/cf. The yields were obtained within 
fiducial and EC cuts. An e^e~ correction was applied to the 
data. No acceptance or efficiency correction was applied to 
either spectrum. 
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FIG. 7: Normalized yield of elastic electron scattering off the 
deuteron at = 0.47 — 0.63 (GeV/c)^: the data are shown 
by the triangles and the simulations are shown by the open 
circles. No acceptance or efficiency corrections were applied. 



the differences obtained for the efficiency were taken as an 
estimate of the systematic uncertainty. These systematic 
uncertainties were averaged over both kinematic vari- 
ables to give a uniform systematic uncertainty to both 
data sets, which was estimated to be 1.7%. The sec- 
ond systematic uncertainty is due to the inabihty of the 
GEANT3-based CLAS simulation package GSIM HH to 
perfectly reproduce the CLAS response to electron tracks 
at different angles and momenta. To estimate this effect, 
we treated the six CLAS sectors as independent spec- 
trometers. The normalized event yield measured in each 
sector was compared separately to the simulations, as 
shown in Fig. |S| The observed differences were compared 
sector-by-sector to remove uncertainties due to the event 
generator model. From this comparison we obtained a 
systematic uncertainty varying from 3 to 6% depending 
mostly on the scattered electron polar angle. The two 
uncertainties were summed in quadrature. 



E. Structure Function F2{x,Q^) 

The measured electron yields N^xp, normalized to the 
integrated luminosity in conjunction with Monte Carlo 
simulations, were used to extract the structure function 
F2 in each kinematic bin. The Monte Carlo events were 
used to simultaneously obtain efficiency, acceptance, bin 
centering and radiative corrections. F2 was determined 
using: 



1 



J 



PT^^Qtot ^Mott 1 + ^Ttn 



l-e 1 



^^^^F,,,(.T, Q2)F,+ ,- [x, Q') , (15) 

where p is the density of liquid D2 in the target, Na is 
the Avogadro constant, Ma is the target molar mass, L is 
the target length, Qtot is the total charge in the Faraday 
Cup (FC) and ^'(a;,(5^) is the efficiency including the 
radiative and bin-centering correction factors: 

^ix,Q^)^-^effix,Q^)-^radix,Q^)^Mnix,Q^) , (16) 

with: 



and *b.;„ = 



M 



(17) 



of the scattered electron and deuteron, respectively; 
and Pd are the momenta of the particles. As one can see 
from Fig. [7| the elastic peak is very well separated from 
the inelastic background. The integrated peak yields 
agree to within the 3% statistical uncertainty. 

There are two systematic uncertainties in the simula- 
tion. The first one is due to the model dependence of 
the reaction cross section used for generating the events. 
We applied a different cross section model for the inelas- 
tic electron-deuteron scattering, taken from Ref. [i^l , and 



^eff is the ratio between the number of reconstructed 
and generated events in the bin. The integral in Eq. 1171 
was taken over the current bin area Ar. Here e is the 
virtual photon polarization parameter: 



l + 2^tan^- 



(18) 



The Mott cross section aMott and the Jacobian J of the 
transformation between dfldE' to dxdQ^ are defined by: 
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OMott 



o? cos^ I 
4£;2 sin^ % 



and J 



xEE' 



(19) 



The structure function F2(x^ Q^) was extracted using the 
fit of the function R{x,Q^) = OL/dT described in Ap- 
pendix^ However, the structure function F2 in the rele- 
vant kinematic range is very insensitive to the value of R. 
For example in typical kinematics e = 0.75 and assum- 
ing SLAC DIS R value of 0.18 the relative uncertainties 
of F2 and R are related: dFa/Fa = OmdR/R. There- 
fore, in this kinematics 20% error on R will generate only 
0.6% error on F2. The overlapped data from two differ- 
ent beam energies were combined using weighted average 
technique. Moreover, we checked that the parameteriza- 
tion used for R is consistent with the difference between 
two cross sections within statistical and systematic er- 
rors. 

Fig. |S1 shows a comparison between the F2 data from 
CLAS and the other world data in a few bins. The 
CLAS data agree very well with all previous measure- 
ments. The values of F2{x,Q^), together with their sta- 
tistical and systematic uncertainties, are tabulated else- 
where 42]. 

In the calculation of the radiative correction factor 
^rad we used the cross section model described in Ap- 
pendix ^ in the following way: 

• the eD elastic radiative tail was calculated accord- 
ing to the "exact" Mo and Tsai formula |3^ : 

• in the quasi-elastic peak region {W'^'' + AW < W < 
1.2 GeV) we applied the correction formula to the 
continuum spectrum given in Ref. j36(|, which is 
based on the peaking approximation and is known 
to be rehable only when E'/E > 0.5. Here W^^ is 
the eD elastic peak position and AW its width; 

• at > 1.2 GeV we applied the exact Mo and Tsai 
formula also to the quasi-elastic tail and a peak- 
ing approximation based formula (referred to as 
the "unfolding procedure") to the inelastic spec- 
trum. For an exact calculation of the quasi-elastic 
tail it was necessary to extract quasi-elastic form- 
factors. To this end we integrated the quasi-elastic 
cross section from the beginning of the peak up to 
W = 1.2 GeV and performed a separation of the 
electric and magnetic form-factors. These two kine- 
matic regions overlap quite well and don't exhibit 
any discontinuity at the point W = 1.2 GeV. This 
assured us that the peaking approximation formal- 
ism is safely applicable to the quasi-elastic tail up 
toW = 1.2 GeV. 

The radiative correction factor "^rad varies strongly in 
the explored kinematic range from 0.7 up to 1.5. Fortu- 
nately, the largest corrections are given by the tails of the 
elastic and quasi-elastic peaks, for which calculations are 



very accurate (see Refs. Js^, 'i^)- The largest system- 
atic uncertainties (see Table are due to the efficiency 
evaluation, and to the photo-electron and radiative cor- 
rections. 

The systematic uncertainties for the efficiency evalu- 
ation, the e^e~ pair correction, and the pion rejection 
correction were described above, and the systematic un- 
certainties arising from the applied CLAS momentum 
correction routines are calculated according to Ref. Q. 
The radiative correction factors in cr^f^ were evaluated 
with two different methods ( 36] , 43] ) and the difference 
was taken as an estimate of the corresponding systematic 
uncertainty. These two methods use different parameter- 
izations of the elastic (|3^ and li^ll . quasi-elastic ( [iol ] 
and ^5) and inelastic (3 and ^4^) cross sections, as 
well as different calculation techniques. The uncertain- 
ties in R given in Appendix were propagated to the 
resulting F2. All systematic uncertainties were summed 
in quadrature to obtain the final systematic uncertainty. 

The statistical and systematic precisions of the ex- 
tracted structure function F2 are strongly dependent on 
the kinematics: the statistical uncertainties vary from 
0.1% up to 30% at the largest Q^, where the event yield 
is very limited, while the average value is about 3%; the 
systematic uncertainties range from 4% up to 14%, with 
the mean value being about 7% (see Table QJ. 



TABLE I: Range and average of systematic uncertainties on 
F2. 



Source of uncertainties 


Variation range 


Average 


]%] 


]%] 


Efficiency evaluation 


3-7 


4.6 


e^e~ pair production correction 


0-3 


0.1 


Pion rejection correction 


0.1-6 


2.4 


Radiative correction 


1.8-3.5 


2.3 


Momentum correction 


0.2-1.2 


0.5 


Uncertainty of i? = ^ 


0.2-0.75 


0.5 


Empty target subtraction 


0.4-0.42 


0.4 


Total 


4-14 


6.6 



F. Moments of the Structure Function F2 

The evaluation of the deuteron structure function mo- 
ments was performed according to the method developed 
m Ref. ^. However, there are two main differences in 
the deuteron data analysis: 

1. the quasi-elastic peak is not as well known as the 
proton elastic form-factors, and moreover cannot 
be easily separated from the inelastic spectrum. 
Hence, in contrast to the free proton target, we ex- 
tract the total moments of the deuteron structure 
function F2 directly, without separating them into 
the elastic and inelastic parts. This emphasizes the 



11 




FIG. 8: The deuteron structure function F2(x,Q^) per nucleon at six different values: from top-left in order = 0.825 
(GeV/c)^ = 1.375 (GeV/c)^ = 1.775 {GeY/cf. = 2.825 (GeV/c)^ = 4.075 (GeV/c)^ and = 5.025 
(GeV/c)^. The triangles represent experimental data obtained in the present analysis with systematic uncertainties indicated 
by the hatched area. The empty circles show data from previous experiments li iS, iS, ili lUi, il^ ili ili il^ il^ • 
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importance of a precise determination of the quasi- 
elastic peak for each together with the inelastic 
spectrum. In particular, the contribution of the 
quasi-elastic peak in the higher moments (n > 2) 
at values in the interval 1-5 (GeV/c)^ is very 
significant, as one can see in Fig.lHl Thanks to the 
CLAS data this problem is well addressed now; 

2. the lack of collider data does not allow one to reach 
very low x values. For the second moment M2 this 
leads to an increase of systematic uncertainties due 
to the low-x extrapolation with respect to the pro- 
ton M2. The contribution of the low x part in the 
higher moments however is negligible. We used two 
models of the deuteron structure function F2 which 
have very different low- a; behaviour to estimate the 
extrapolated part of the second moment M2. The 
difference between the two estimates was taken as 
an evaluation of the corresponding systematic un- 
certainty. A comparison of the extrapolation sys- 
tematic uncertainties and other uncertainties in M2 
is shown in Fig. IIUI 

We combined the structure functions F2 obtained from 
the CLAS data and the other world data on the struc- 
ture function F2, alongwith the inclusive cross section 
data from RefsJl, HBllllSlillllllllll 111 13 
ElllllllllSIIJlsee Fig.P). The data from Ref. Hj, 
recently reanalyzed in Ref. .45] with the inclusion of ra- 
diative and bin centering corrections, are not used in the 
present analysis due to large statistical uncertainties and 
unknown systematic uncertainties. The Q^-range of the 
CLAS data, from 0.4 to 5.95 (GeV/c)^, was divided into 
bins of width AQ2 = o.05 (GeV/c)^. Within each 
bin the world data were shifted to the central bin value 
Ql, using the fit of Fi^{x,Q^) from Appendix. El The 
integrals of the data over x were performed numerically 
using the standard trapezoidal method TRAPER [4y|. 
As an example. Fig. shows the integrands of the first 
four moments as a function of x at fixed Q^. The sig- 
nificance of the large x region for various moments can 
clearly be seen. 

As in Ref. ||, the world data at above 6 (GeV/c)^ 
were analyzed in the same way as described above, but 
with a different bin size. The bin size was chosen for 
the data to provide sufficient x-coverage for most of the 
bins (AQ2/Q2 = 5 %) The results together with 
their statistical and systematic uncertainties are shown 
in Fig.lTTland reported in Table ITTl 

The systematic uncertainty consists of experimental 

uncertainties in the data given in Refs. [illki,ia,il 

ElllllllllllllllllillllllSilland uncertainties in 
the evaluation procedure. To estimate the first type of 
uncertainty we had to account for the inclusion of many 
data sets measured in different laboratories with differ- 
ent detectors. In the present analysis we assume that the 
different experiments are independent and therefore only 
the systematic uncertainties within a given data set are 
correlated. 
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X 

FIG. 9: Integrands of the Nachtmann moments at — 0.825 
(GeV/c)^: circles represent the integrand of the M2; squares 
show the integrand of the M4; triangles show the integrand 
of the Afe; crosses show the integrand of the Mg. 
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FIG. 10: Uncertainties of the Nachtmann moment M2 in per- 
centage. The lower cross-hatched area represents statistical 
uncertainties. The left-hatched area represents the systematic 
uncertainties. The right-hatched area represents the low- a; ex- 
trapolation uncertainty. 



IV. SEPARATION OF LEADING AND HIGHER 
TWISTS 

In order to separate the leading and higher twists in 
the measured moments we used the method developed 
in Refs. 4, 24, 47] . This method is essentially based 
on a general form of the OPE for the structure function 
moments, where the leading twist Q^-evolution is calcu- 
lated in pQCD and the deviation from this behaviour is 
assigned to the higher twist contribution. Therefore the 
measured Nachtmann n-th moment was parameterized 
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where A„ = /3oas(Q^)ln(n)/47r and: 
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FIG. 11: The Nachtman moments extracted from the world 
data, including the new CLAS results. Uncertainties are sta- 
tistical only. 



as: 



(20) 



where rin{Q^) is the leading twist moment and HTn(Q^) 
is the higher twist contribution. The leading twist term 
was calculated including soft giuon re-summation (SGR) 
corrections which go beyond the fixed-order next-to- 
leading approximation and are essential for a reliable ex- 
traction of higher twists [2^. The observed decoup ling 
of the singlet quark and gluon densities at large x |43| 
allows us to consider only the non-singlet (NS) evolution 
for n > 4 and therefore to reduce the number of lead- 
ing twist parameters. In this approximation the leading 
twist moment rjniQ^) for n > A can be written as follows: 



-R: 



NS 



(21) 



where the quantities 7^"^, c\^jg^^ and R^'^ can be ob- 
tained from Ref. [H, a,(M|) = 0.118 @ and the ref- 
erence scale fM^ = 10 (GeV/c)^. In Eq. |2] the func- 
tion Gn{Q^) is the key quantity of the soft gluon re- 
summation. At Next-to-Lcading-Log it reads as: 



G„(Q2) = ln(n)Gi(A„) + G2(A„) + 0[a^ln'=-i (n)] ,(22) 



Gi(A) 
G2(A) 



Cf-ti- [A +(1-A)ln(l^ A)] . 
PoA 

47B + 3 
-Cf — -s ln(l - A) 

Po 



-GF^[A-hln(l-A)] 



-Gj 



8K 



A + ln(l-A) + iln2(l-A) 



(23) 



with Cf = {N^ ~ l)/(2A^c), k = Ar^(67/18 - - 
5iV//9, /3o = 11 - 2iV//3, and Nf being the number of 
active flavors. Note that the function G2(A) is diver- 
gent for A — > 1. This means that at large n (i.e. large 
x) the soft gluon re-summation cannot be extended to 
arbitrarily low values of Q^. Therefore, for a safe use 
of present SGR techniques we work far from the above- 
mentioned divergences by limiting our analysis of low- 
order moments (n < 8) to > qj - 1 (GeV/c)2. 

Since a complete calculation of the higher twist anoma- 
lous dimensions is not yet available, we use the same phe- 
nomenological ansatz already adopted in Refs. ■4 i24i47j |. 
In this approach the higher twist contribution is given 
by M- 



.(4) 



.(6) 





7'-" 2 








7<'' 4 







(24) 



where the logarithmic pQCD evolution of the 
twist-T contribution is accounted for by the term 

[as(Q^)/as(A*^)]^"^' • This term corresponds to the 
Wilson coefficient ii^„T-(/ir, /i/, Q^) in Eq. |21 with an 

effective anomalous dimension \ The parameter an ^ 
represents the overall strength of the twist-r term at the 
renormalization scale = /i^ and it is proportional to 
the matrix element 0„r (/^) in Eq- [21 The presence of two 
distinct higher twist terms for n > 4 is motivated by the 
Q^-behaviour of the total higher twist contribution in 
the moments. This was obtained by a direct subtraction 
of the leading twist term fitted to the large part of the 
plot from the measured moments. Existence of maxima 
and moreover of the sign turn-over (see Fig. I12|) in the 
total higher twist contribution cannot be described by a 
single twist term within the pQCD-inspired model from 
Eq. 1^ Therefore, the presence of at least two higher 
twist terms is necessary for a successful description of 
experimental data. We checked that the variation of the 
total higher twist contribution after inclusion of twist-8 
and twist- 10 terms is smaller than the quoted systematic 
uncertainties. 

The n-th moment (see Eqs. EOl |^ and |^ for n > 4 

has five unknown parameters: the twist-2 parameter 
An and the higher twist parameters a'n' , ^If' , ai^"* , 7^^^ . 
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All five unknown parameters were simultaneously deter- 
mined fi'om a x^-minimization procedure in the range 
between 1 and 100 (GeV/c)'^. In this procedure only 
the statistical uncertainties of the experimental moments 
were taken into account. The uncertainties of the various 
twist parameters were then obtained by adding the sys- 
tematic uncertainties to the experimental moments and 
by repeating the twist extraction procedure. 

For n = 2, ?72(Q^) is given by the sum of the non-singlet 
and singlet terms, which yield two unknown parameters 
associated with the leading twist. These parameters are 
the values of the gluon and non-singlet quark moments 
at the reference scale — However, due to the 
vanishing contribution of the higher twists in (see 
Fig. I12() one can reduce the number of parameters in 
HTn{Q^) by limiting the expansion to the twist-4 term 
only. 

The parameter values obtained at the reference scale 
^2 = iQ (GeV/c)2 are reported in Table IIIII where it 
can be seen that the leading twist is determined with 
an uncertainty of a few percent, while the precision of 
the extracted higher twists decreases with n, reaching an 
overaU 20-30% for n = 8, thanks to the CLAS data at 
large x. Note that the leading twist is directly extracted 
from the data, which means that no specific functional 
shape of the parton distributions is assumed in our anal- 
ysis. 

Our results for each twist term are reported in Fig. 1121 
for 71 > 2, while the ratio of the total higher twist con- 
tribution to the leading twist is shown in Fig. 1131 In 
addition, the extracted leading twist contribution is re- 
ported in Table HVl 



V. CONCLUSIONS 

We extracted the deuteron F2 structure function in 
a continuous two-dimensional range of and x from 
the inclusive cross sections measured with CLAS. The 
extracted structure functions are in good agreement 
with previous measurements in overlapping regions, con- 
tribute many additional kinematic points and sometimes 
improve the precision where the world data exist. Us- 
ing these data, together with the previously available 
world data set, we evaluated for the first time the Nacht- 
mann moments M2{Q^,x), Mi{Q^ ,x), Mq{Q^,x) and 
MsiQ^,x) in the range 0.5 - 100 (GeV/c)2. Pre- 
viously, the experimental information on the deuteron 
structure function moments was missing in the low to 
medium domain due to scarce data at a; ^ 1. More- 
over, fixed bins of the data render the Q'^-evolution 
of the extracted moments model- independent. The mea- 
sured moments have been analyzed in terms of a twist 
expansion in order to extract both the leading and the 
higher twists simultaneously. By calculating the Q^- 
evolution of the leading twist at NLO and including 
Qfs-higher order corrections through the soft gluon re- 
summation, we extracted values of the reduced matrix 



elements 0„2 from Eq.|21 Higher twists have been treated 
phenomenologically within the pQCD-inspired approach 
of Eq. 121 by introducing effective anomalous dimensions. 
The interval of the analysis was quite large, ranging 
from 1 to 100 (GeV/c)^, allowing us to determine the 
total contribution of higher twists to the best accuracy 
possible. The variation of the total higher twist con- 
tribution due to inclusion of twist-8 and twist- 10 terms 
is lower than the quoted systematic uncertainties. The 
leading twist is determined with a few percent uncer- 
tainty, while the precision of the higher twists decreases 
with n reaching an overall 20-30% for n =8, thanks to 
the remarkable quality of the experimental moments. 

The main results of our twist analysis can be summa- 
rized as follows: 

• the extracted leading term yields the dominant con- 
tribution in the entire Q^-range of the present anal- 
ysis for all four moments. This leads to the conclu- 
sion that despite the nuclear effects, a pQCD-based 
description of the deuteron structure, including the 
effects of soft-gluon re-summation, is surprisingly 
applicable also at low Q^. The corrections to the 
leading twist are significant but not crucial; 

• the Q^-behaviour of the data indicates the presence 
of the higher twist contribution at <5 (GeV/c)^, 
positive at large and negative at ^1- 
2 (GeV/c)^; the change of sign requires in Ea. 1241 
at least two higher twist terms with opposite signs. 
As already noted in Refs. 0, 01 such a can- 
cellation makes the total higher twist contribution 
smaller than its individual terms, which exceed the 
leading twist. This partial cancellation is a man- 
ifestation of the duality phenomena in the pQCD 
representation |50j |. It leads to the prevailing DIS- 
inspired picture of virtual photon- nucleon collisions 
also at low Q^. The same mutual cancellation of 
higher twist terms was observed in the proton struc- 
ture function moments in Refs. [l^lsij: 

• the total higher twist contribution is significant at 

« few (GeV/c)^ and large x. This can be seen 
by comparing the higher twist contribution to Mg, 
which is more heavily weighted in x, to M2. For 

>6 (GeV/c)2 the higher twist contribution does 
not exceed ~20% of the leading twist for all four 
moments. 

Therefore, we have demonstrated that despite nuclear 
effects in the deuteron, a pQCD-based analysis of the 
deuteron structure function moments is sensible, so that 
a precise determination of the leading and higher twists is 
possible with the new CLAS data. The extracted values 
of the reduced matrix elements 0„2 still contain some 
contribution of the nuclear off-shell and Fermi motion 
effects, which should be taken care of before a compari- 
son to Lattice QCD simulations of the nucleon is made. 
However, most of the nuclear effects, in particular FSI, 
should be absorbed in the effective higher twist terms 
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TABLE II: The Nachtmann moments for n = 2, 4, 6 and 8 evaluated in the interval 0.05 < < 100 (GeV/c)^. The moments 
are labeled with an asterisk when the contribution to the integral by the experimental data is between 50% and 70%. All the 
others were evaluated with more than 70% data coverage. The data are reported together with the statistical and systematic 
uncertainties, the third uncertainty for n = 2 is due to low-a; extrapolation. 



[(GoV/c)"] 


M2(Q^)xl0-' 


M4(Q^)xl0-^ 


Af6(Q^)xlO-^ 


M8(Q^)xlO-'' 


0.475 


2.133 ± 0.014 ± 0.063 ± 0.004 


3.673 ± 0.035 ± 0.139 


0.911 ± 0.010 ± 0.037 


2.421 ± 0.030 ± 0.101 


0.525 


2.168 ± 0.005 ± 0.093 ± 0.003 


3.868 ± 0.017 ± 0.171 


1.011 ± 0.005 ± 0.046 


2.852 ± 0.016 ± 0.129 


0.575 


2.118 ± 0.004 ± 0.082 ± 0.001 


3.776 ± 0.011 ± 0.168 


1.026 ± 0.004 ± 0.047 


3.042 ± 0.012 ± 0.140 


0.625 


2.084 ± 0.005 ± 0.073 ± 0.003 


3.757 ± 0.008 ± 0.162 


1.067 ± 0.003 ± 0.048 


3.328 ± 0.010 ± 0.154 


0.675 


2.064 ± 0.006 ± 0.077 ± 0.003 


3.757 ± 0.009 ± 0.179 


1.101 ± 0.003 ± 0.053 


3.572 ± 0.013 ± 0.172 


0.725 


2.027 ± 0.011 ± 0.067 ± 0.002 


3.632 ± 0.010 ± 0.168 


1.086 ± 0.004 ± 0.053 


3.641 ± 0.016 ± 0.178 


0.775 


1.970 ± 0.008 ± 0.066 ± 0.003 


3.515 ± 0.012 ± 0.155 


1.080 ± 0.005 ± 0.051 


3.752 ± 0.022 ± 0.182 


0.825 


1.955 ± 0.010 ± 0.073 ± 0.001 


3.460 ± 0.012 ± 0.174 


1.079 ± 0.005 ± 0.055 


3.840 ± 0.024 ± 0.195 


0.875 


i.9()8 ± 0.005 ± O.OGl ± 0.001 


:!. 171 ± 0.007 ± 0.175 


i.l08 ± 0.003 ± 0.058 


l.()5G ^ 0.015 ± 0.213 


0.925 


1.926 ± 0.010 ± 0.062 ± 0.003 


3.357 ± 0.014 ± 0.166 


1.082 ± 0.007 ± 0.057 


4.042 ± 0.035 ± 0.213 


0.975 


1.916 ± 0.011 ± 0.063 ± 0.002 


3.352 ± 0.013 ± 0.162 


1.106 ± 0.006 ± 0.058 


4.238 ± 0.032 ± 0.225 


1.025 




3.280 ± 0.013 ± 0.169 


1.099 ± 0.006 ± 0.060 


4.320 ± 0.034 ± 0.239 


1.075 


1.873 ± 0.009 ± 0.066 ± 0.005 


3.159 ± 0.014 ± 0.164 


1.059 ± 0.007 ± 0.059 


4.215 ± 0.039 ± 0.240 


1.125 


1.837 ± 0.006 ± 0.059 ± 0.006 


3.071 ± 0.007 ± 0.147 


1.036 ± 0.003 ± 0.056 


4.190 ± 0.018 ± 0.235 


1.175 




3.012 ± 0.008 ± 0.132 


1.022 ± 0.003 ± 0.052 


4.201 ± 0.017 ± 0.228 


1.225 




2.934 ± 0.016 ± 0.148 


0.997 ± 0.007 ± 0.055 


4.136 ± 0.037 ± 0.234 


1.275 




2.842 ± 0.013 ± 0.147 


0.965 ± 0.006 ± 0.053 


4.035 ± 0.032 ± 0.230 


1.325 




2.809 ± 0.012 ± 0.154 


0.961 ± 0.005 ± 0.055 


4.080 ± 0.027 ± 0.238 


1.375 


1.807 ± 0.008 ± 0.055 ± 0.009 


2.789 ± 0.014 ± 0.138 


0.951 ± 0.008 ± 0.053 


4.064 ± 0.046 ± 0.238 


1.425 




2.751 ± 0.011 ± 0.130 


0.944 ± 0.006 ± 0.051 


4.086 ± 0.033 ± 0.235 


1.475 




2.681 ± 0.016 ± 0.119 


0.926 ± 0.007 ± 0.047 


4.063 ± 0.041 ± 0.221 
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2.582 ± 0.011 ± 0.140 


0.894 ± 0.005 ± 0.052 


3.992 ± 0.025 ± 0.236 


1.675 


1.745 ± 0.011 ± 0.069 ± 0.019 


2.575 ± 0.013 ± 0.135 


0.891 ± 0.007 ± 0.051 


4.012 ± 0.042 ± 0.239 


1.725 


1.733 ± 0.006 ± 0.046 ± 0.011 


2.527 ± 0.008 ± 0.115 


0.872 ± 0.005 ± 0.046 


3.942 ± 0.033 ± 0.225 


1.775 


1.721 ± 0.010 ± 0.055 ± 0.019 


2.468 ± 0.009 ± 0.109 


0.845 ± 0.006 ± 0.041 


3.816 ± 0.035 ± 0.200 


1.825 


1.709 ± 0.010 ± 0.054 ± 0.021 


2.413 ± 0.010 ± 0.111 


0.829 ± 0.007 ± 0.041 


3.789 ± 0.048 ± 0.193 


1.875 




2.402 ± 0.013 ± 0.117 


0.824 ± 0.006 ± 0.045 


3.774 ± 0.031 ± 0.212 


1.925 




2.349 ± 0.018 ± 0.112 


0.794 ± 0.008 ± 0.041 


3.610 ± 0.047 ± 0.200 


1.975 


1.697 ± 0.014 ± 0.049 ± 0.023 


2.315 ± 0.012 ± 0.105 


0.777 ± 0.008 ± 0.039 


3.519 ± 0.060 ± 0.186 


2.025 


1.698 ± 0.007 ± 0.052 ± 0.011 


2.280 ± 0.005 ± 0.106 


0.762 ± 0.003 ± 0.038 


3.456 ± 0.018 ± 0.178 


2.075 




2.246 ± 0.011 ± 0.104 


0.747 ± 0.004 ± 0.037 


3.395 ± 0.023 ± 0.171 


2.125 


1.696 ± 0.016 ± 0.061 ± 0.023 


2.244 ± 0.010 ± 0.110 


0.744 ± 0.007 ± 0.038 


3.383 ± 0.054 ± 0.172 


2.175 


1.684 ± 0.012 ± 0.066 ± 0.023 


2.211 ± 0.010 ± 0.109 


0.731 ± 0.007 ± 0.037 


3.324 ± 0.052 ± 0.172 


2.225 


1.687 ± 0.009 ± 0.044 ± 0.019 


2.201 ± 0.004 ± 0.103 


0.728 ± 0.002 ± 0.036 


3.324 ± 0.016 ± 0.168 


2.275 




2.168 ± 0.019 ± 0.099 


0.716 ± 0.004 ± 0.035 


3.274 ± 0.023 ± 0.162 


2.325 


1.683 ± 0.011 ± 0.055 ± 0.025 


2.148 ± 0.006 ± 0.094 


0.700 ± 0.004 ± 0.032 


3.179 ± 0.029 ± 0.149 


2.375 


1.669 ± 0.010 ± 0.041 ± 0.032 


2.104 ± 0.005 ± 0.085 


0.685 ± 0.003 ± 0.030 


3.121 ± 0.027 ± 0.143 


2.125 




2.003 ± O.OOG ± 0.088 


O.OOG ± 0.003 ± 0.028 


3.01G I 0.021 ± 0.130 


2.475 




2.054 ± 0.008 ± 0.079 


0.661 ± 0.003 ± 0.027 


2.991 ± 0.020 ± 0.125 


2.525 


1.639 ± 0.026 ± 0.059 ± 0.029 


2.034 ± 0.005 ± 0.089 


0.653 ± 0.002 ± 0.028 


2.956 ± 0.019 ± 0.121 


2.575 


1.644 ± 0.013 ± 0.051 ± 0.027 


2.032 ± 0.008 ± 0.091 


0.656 ± 0.005 ± 0.032 


2.990 ± 0.042 ± 0.144 


2.625 




2.026 ± 0.009 ± 0.101 


0.652 ± 0.004 ± 0.034 


2.974 ± 0.032 ± 0.153 


2.675 




1.986 ± 0.008 ± 0.090 


0.637 ± 0.003 ± 0.032 


2.914 ± 0.021 ± 0.154 


2.725 


1.652 ± 0.009 ± 0.040 ± 0.022 


1.990 ± 0.004 ± 0.093 


0.630 ± 0.002 ± 0.032 


2.865 ± 0.017 ± 0.152 


2.775 


1.620 ± 0.019 ± 0.042 ± 0.016 


1.961 ± 0.015 ± 0.089 


0.623 ± 0.008 ± 0.032 


2.819 ± 0.062 ± 0.151 


2.825 




1.939 ± 0.029 ± 0.092 


0.613 ± 0.007 ± 0.032 


2.780 ± 0.038 ± 0.149 


2.875 




1.912 ± 0.015 ± 0.103 


0.600 ± 0.008 ± 0.031 


2.707 ± 0.060 ± 0.144 


2.925 




1.891 ± 0.005 ± 0.090 


0.591 ± 0.003 ± 0.030 


2.668 ± 0.020 ± 0.139 


2.975 




1.896 ± 0.016 ± 0.096 


0.588 ± 0.007 ± 0.032 


2.641 ± 0.057 ± 0.142 


3.025 


i.()21 ± 0.02i ± 0.015 ± 0.031 


i.873 ± 0.011 ± 0.090 


0.581 ± ().00() ± 0.032 


2.G07 ^ 0.057 ± 0.144 


3.075 




1.853 ± 0.039 ± 0.090 


0.571 ± 0.007 ± 0.032 


2.557 ± 0.038 ± 0.145 


3.125 




1.851 ± 0.006 ± 0.100 


0.571 ± 0.003 ± 0.032 


2.566 ± 0.020 ± 0.147 


3.175 




1.852 ± 0.031 ± 0.110 


0.567 ± 0.009 ± 0.034 


2.536 ± 0.059 ± 0.149 


3.225 




1.823 ± 0.007 ± 0.097 


0.557 ± 0.003 ± 0.032 


2.482 ± 0.021 ± 0.146 


:!.275 




1.818 ± 0.018 ± 0.098 


0.551 ± 0.009 ± 0.():i:5 


2.129 = 0.()()7 ± 0.114 
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O.U ( 








1 7S1 4- 01 ft 4- 1 Pil 
i./oi It U.Ulo It U.lOi 


1 9Ci 

O.iZO 




1 /IQQ -L O 01 7 -1- O OQ9 

i.4yy ^ u.ui( ^ U.UOZ 




1 7Kfi 4- O^iQ 4- 1/1/^ 
l./OD ^ u.uoy i U.i40 


1 7t^ 
O.i (0 








1 7c;q j_ 01/^ -1- loci 
l./Oo dl U.UIO it U.lzo 


o.zzo 








1 77 A 4- O'^O 4- 1 '^/l 
i.//4 zt U.UOU zt U.104 


o.z / 






fl/llft 4- 0010 -1- 0097 
U.4io zt U.UiU zt U.UZ( 


1 7Q7 4- OQft 4- 1 Q(^ 
l./o/ zt U.Uoo zt U.loO 


o.ozo 






0/110 -I- 001/1 -1- 009/1 
U.4iy zt U.Ui4 zt U.UZ4 


1 7QQ 4- O^iO 4-0 1 QQ 

i./oy zt U.uoy zt u.ioo 


0.0 / 1) 






n lie; -u o oof. j- o ooi 

U.4il) zt U.UUo zt U.UZ4 


1 710 4- 09 1 -k 1 9 Pi 
l./lU zt U.UZI zt U.lZO 










1 ~(U: — /I (VJ(l -U n 1 10 
l./Uu ziz U.U'jU zt U.14y 


0.4 ( 


1.040 ^ U.UiO ^ U.UoZ ^ U.Uoo 






1 70/1 -1- O^i/I -1- 1 Pi /I 
1./U4 ^ U.U04 i U.i04 


O.OZO 








1 fi/17 J_ OQQ -1- lAA 
1.04/ ± U.Uoo ± U.144 


o.ozo 






O Af\A 4- O OOQ -1- O 090 
U.4U4 ^ u.uuo ^ U.UZU 


1 (^/1Q 4- 000 4- 1 1 Q 
i.D40 zt U.UZU zt U.iiO 


o.yzo 








1 i^OQ 4- 00/^ 4- 00ft 

i.Duo zt U.UZO zt U.uyo 


o.yoo 




1 /IQ^i 4- O 09Q -1- 097 
1.400 It U.UZO zt U.UZi 


n Q7/1 -U 00ft -1- 01 1 
U.0/4 zt U.UUo zt U.Uii 


1 /170 4- OQO 4- O^iO 
i.4/Z zt U.UOZ zt U.UOZ 


o.yio 


i.ozi ^ u.uuy ^ u.ui4 ^ U.UOO 


1 Af\A -1- O 01 1 -1- O 097 

i.4U4 ^ U.Uii ^ u.uzr 


Qf^l -1- O 00/1 -1- O 01 O 
U.OOi ^ U.UU4 ^ U.UiU 


1 QSQ -1- 000 -1- 0/1/1 
l.Ooy ^ U.UZU ^ U.U44 


7 9(^7 

r .ZD ( 




1 Q7f; -1- O m 9 -1- O OQfl 

i.o/0 ± U.Uiz ± U.UoO 


o qc;q -1- o 00/1 -1- o 010 
U.ooo ± U.UU4 ± U.Ulz 


1 QAQ J_ 000 -1- O/^c; 

1.000 dl U.Uzz ± U.UOO 


7 fi'^n 

/ .DOU 






n Q/1Q 4- 00/1 -1- 01 

U.040 zt U.UU4 zt u.uiy 


1 QOS 4- 000 4- lOPi 
l.OUo zt U.UZZ zt U.lZO 


o.UZl 






n QQ/^ -1- 009 4-0 019 
u.ooo zt U.UUZ zt U.UiZ 


1 07/1 4- 010 4- OPiO 
i.Z/4 zt U.UiU zt U.UOZ 


ft S/17 


l.OUo ^ U.UiO ^ U.UiO ^ U.U4i 


1 QOC; -1- 01 1 -1- 097 
l.OZO zt U.Ull zt U.UZi 


n Q90 4- OOQ 4- 010 

u.ozy zt u.uuo zt U.UiU 


1 OlPi 4- 0017 4- 00/1/1 
i.ZiO zt U.Ui/ zt U.U44 


Q 77Ci 

y. M 




1 OQ1 _4_ n OO^i -U n OQS 
l.Zol It U.UUO It U.UOo 


Q1 Q 4- O 001 -1- O 01 O 

u.oio di u.uui ^ U.UiU 


1 1/lR 4- OOPi -1- 0/1 1 
1.140 dl u.uuo ^ U.U4i 


1 n 9*^7 
iU.zu r 








1 1 J_ 01 -1- OQ1 
l.lOu ± U.UIU ± U.Uol 


10.762 






0.306 lb 0.002 lb 0.007 


1.115 lb 0.011 lb 0.029 


11.344 


1.500 ± 0.014 ± 0.010 ± 0.046 


1.242 ± 0.015 ± 0.021 


0.299 lb 0.005 lb 0.008 


1.084 lb 0.035 ± 0.033 


12.580 








1.049 zb 0.013 zb 0.022 


13.238 


1.503 ± 0.025 ± 0.023 ± 0.049 


1.205 ± 0.011 ± 0.019 


0.287 ± 0.003 zb 0.006 


1.001 zb 0.010 zb 0.026 


14.689 


1.520 ± 0.012 ± 0.011 ± 0.052 


1.216 ± 0.013 ± 0.027 


0.285 ± 0.003 zb 0.010 


0.992 zb 0.013 zb 0.048 


17.108 


1.478 ± 0.011 ± 0.017 ± 0.052 


1.112 ± 0.011 ± 0.038 


0.256 ± 0.003 ± 0.011 


0.900 ± 0.015 lb 0.039 
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[{GeV/cY] 


M2(Q^)xlO~' 


M4(Q^)xlO"" 


M6(0^)xlO"" 


Af8(Q")xlO-^ 


19.072 




1.107 ± 0.010 ± 0.021 


0.257 ± 0.003 ± 0.007 


0.911 ± 0.014 ± 0.028 


20.108 


1.467 ± 0.010 ± 0.007 ± 0.058 


1.094 ± 0.014 ± 0.019 


0.250 ± 0.004 ± 0.008 


0.874 ± 0.016 ± 0.039 


21.097 


*1.486 ± 0.010 ± 0.011 ± 0.049 


1.101 ± 0.012 ± 0.032 






24.259 


*1.470 ± 0.012 ± 0.014 ± 0.051 


1.052 ± 0.010 ± 0.040 






26.680 


*1.469 ± 0.013 ± 0.007 ± 0.059 


1.056 ± 0.010 ± 0.021 






32.500 


*1.479 ± 0.013 ± 0.014 ± 0.055 








34.932 


*1.445 ± 0.113 ± 0.007 ± 0.064 








36.750 


*1.482 ± 0.012 ± 0.018 ± 0.057 


1.014 ± 0.010 ± 0.016 


0.223 ± 0.003 ± 0.004 


0.767 ± 0.014 ± 0.014 


43.970 


*1.478 ± 0.026 ± 0.010 ± 0.057 








47.440 


*1.477 ± 0.127 ± 0.007 ± 0.051 


1.016 ± 0.029 ± 0.018 






64.270 




9.450 ± 0.042 ± 0.014 


0.206 ± 0.007 ± 0.003 


0.699 ± 0.022 ± 0.010 


75.000 




*9.443 ± 0.017 ± 0.024 


*0.206 ± 0.006 ± 0.008 




86.000 




*9.183 ± 0.018 ± 0.019 


*0.196 ± 0.007 ± 0.007 




97.690 




*9.232 ± 0.019 ± 0.010 


*0.199 ± 0.005 ± 0.003 


*0.672 ± 0.023 ± 0.008 



TABLE III: Extracted parameters of the twist expansion at the reference scale /i^ — 10 (GeV/c)^. The first uncertaint y is 
the systematic one described in text, while the second uncertainty has a statistical origin and is obtained from a MINOS |4q| 
minimization procedure. The contribution of twists-6 to A/2 was too small to be extracted by the present procedure. 





M2 


M4 


Me 


Ms 




0.152±0.02±0.03 


(1.215±0.03±0.01)xl0"^ 


(2.95±0.1±0.04)xl0~^ 


(1.05±0.02±0.02)xl0"'' 




(4±2±22)xl0"'' 


(7.4±2±2.5)xl0"^ 


(2.7±0.5±0.04)xl0"^ 


(1.7±0.8±0.04)xl0"^ 




3.4±0.2±5.2 


3±0.5±1 


5.9±0.3±0.02 


6.4±3.5±0.04 






(-1.5±0.2±0.3)xl0-^ 


(-9.2±1.4±0.15)xl0~^ 


(-6.6±2±0.16)xl0~^ 


^(6) 




1.9±0.4±0.8 


4.3±0.3±0.02 


4.7±1.4±0.04 



due to their scale difference. An estimate of the leading 
twist nuclear corrections would allow extraction of the 
non-singlet part of the nucleon structure function mo- 
ments, which can be directly tested in the Lattice QCD 
simulations. 
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APPENDIX A: MODEL OF INCLUSIVE 
ELECTRON SCATTERING CROSS SECTION 
OFF THE DEUTERON AND THE DEUTERON 
STRUCTURE FUNCTION F2 

In order to extract efficiency, calculate radiative cor- 
rections and evaluate moments of the structure function 
F2, it is essential to have a realistic model of the reac- 
tion cross section. Thanks to many previous experiments, 



comprehensive knowledge about electron-deuteron inclu- 
sive scattering is available. We based our model on these 
previous results. The model consists of three main ele- 
ments: 

• the elastic peak cross section was calculated using 
the deuteron elastic form-factors from Ref. ,34, .52j ; 

• the quasi-elastic peak cross section was obtained 
within a model of the nuclear structure of the 
deuteron '40'| using elementary form-factors of the 
proton and neutron from Ref. jH^ . This model is 
based on the De Forest 54] ccl prescription for the 
off-shell nucleon and includes various sophisticated 
treatments of the final state interactions. Specifi- 
cally designed for calculations of the quasi-elastic 
cross section for light and complex nuclei, it repro- 
duces existing data very well (see Fig. 1141) ; 

• the inelastic cross section in the CLAS domain 
(ly^ < 4.3 GeV^) was taken from the fit to the 
recent Hall C data P|; 

• in DIS we have chosen the parameterization from 
Ref. , which describes particularly well the low- 
X behaviour of the structure function. 

The elastic peak is very small in our range and hence 
it is only relevant for the radiative correction calculations. 
The inelastic cross section model, which fits the data from 
Hall C very well, was obtained in the same kinematic 
domain. The quasi-elastic cross section calculations are 
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TABLE IV: The extracted leading twist contribution rjn{Q^) (see Eg. I2H shown in Fig. 1121 reported with systematic uncer- 
tainties. 



l(GeV/c)^J 


r?2(Q")xlO-^ 


J74(Q")xlO-" 


r?6(Q")xlO-^ 


r?8(Q")xlO-^ 


1.025 


1.84 ± 0.07 


2.56 ± 0.06 


11.5 ± 0.4 


8.63 ± 0.07 


1.075 


1.83 ± 0.07 


2.49 ± 0.05 


10.6 ± 0.4 


7.28 ± 0.06 


1.125 


1.82 ± 0.07 


2.42 ± 0.05 


9.88 ± 0.3 


6.34 ± 0.05 


1.175 


1.81 ± 0.07 


2.36 ± 0.05 


9.29 ± 0.3 


5.64 ± 0.05 


1.225 


1.80 ± 0.07 


2.31 ± 0.05 


8.80 ± 0.3 


5.11 ± 0.04 


1.275 


1.79 ± 0.07 


2.26 ± 0.05 


8.38 ± 0.3 


4.68 ± 0.04 


1.325 


1.78 ± 0.06 


2.21 ± 0.05 


8.02 ± 0.3 


4.34 ± 0.04 


1.375 


1.77 ± 0.06 


2.17 ± 0.05 


7.70 ± 0.3 


4.05 ± 0.03 


1.425 


1.77 ± 0.06 


2.13 ± 0.05 


7.43 ± 0.3 


3.82 ± 0.03 


1.475 


1.76 ± 0.06 


2.10 ± 0.05 


7.18 ± 0.3 


3.61 ± 0.03 


1.525 


1.75 ± 0.06 


2.07 ± 0.05 


6.96 ± 0.2 


3.43 ± 0.03 


1.575 


1.74 ± 0.06 


2.04 ± 0.04 


6.76 ± 0.2 


3.28 ± 0.03 


1.625 


1.74 ± 0.06 


2.01 ± 0.04 


6.58 ± 0.2 


3.14 ± 0.03 


1.675 


1.73 ± 0.06 


1.98 ± 0.04 


6.41 ± 0.2 


3.02 ± 0.02 


1.725 


1.72 ± 0.06 


1.95 ± 0.04 


6.26 ± 0.2 


2.91 ± 0.02 


1.775 


1.72 ± 0.06 


1.93 ± 0.04 


6.12 ± 0.2 


2.81 ± 0.02 


1.825 


1.71 ± 0.05 


1.91 ± 0.04 


6.00 ± 0.2 


2.73 ± 0.02 


1.875 


1.70 ± 0.05 


1.89 ± 0.04 


5.88 ± 0.2 


2.65 ± 0.02 


1.925 


1.70 ± 0.05 


1.87 ± 0.04 


5.76 ± 0.2 


2.57 ± 0.02 


1.975 


1.69 ± 0.05 


1.85 ± 0.04 


5.66 ± 0.2 


2.50 ± 0.02 


2.025 


1.69 ± 0.05 


1.83 ± 0.04 


5.56 ± 0.2 


2.44 ± 0.02 


2.075 


1.68 ± 0.05 


1.81 ± 0.04 


5.47 ± 0.2 


2.38 ± 0.02 


2.125 


1.68 ± 0.05 


1.80 ± 0.04 


5.39 ± 0.2 


2.33 ± 0.02 


2.175 


1.67 ± 0.05 


1.78 ± 0.04 


5.31 ± 0.2 


2.28 ± 0.02 


2.225 


1.67 ± 0.05 


1.77 ± 0.04 


5.23 ± 0.2 


2.23 ± 0.02 


2.275 


1.66 ± 0.05 


1.75 ± 0.04 


5.16 ± 0.2 


2.19 ± 0.02 


2.325 


1.66 ± 0.05 


1.74 ± 0.04 


5.11 ± 0.2 


2.16 ± 0.02 


2.375 


1.66 ± 0.05 


1.73 ± 0.04 


5.05 ± 0.2 


2.12 ± 0.02 


2.425 


1.66 ± 0.05 


1.72 ± 0.04 


5.00 ± 0.2 


2.09 ± 0.02 


2.475 


1.65 ± 0.05 


1.71 ± 0.04 


4.94 ± 0.2 


2.06 ± 0.02 


2.525 


1.65 ± 0.05 


1.70 ± 0.04 


4.89 ± 0.2 


2.04 ± 0.02 


2.575 


1.65 ± 0.05 


1.69 ± 0.04 


4.85 ± 0.2 


2.01 ± 0.02 


2.625 


1.64 ± 0.05 


1.68 ± 0.04 


4.80 ± 0.2 


1.98 ± 0.02 


2.675 


1.64 ± 0.04 


1.67 ± 0.04 


4.76 ± 0.2 


1.96 ± 0.02 


2.725 


1.64 ± 0.04 


1.66 ± 0.04 


4.71 ± 0.2 


1.93 ± 0.02 


2.775 


1.64 ± 0.04 


1.65 ± 0.04 


4.67 ± 0.2 


1.91 ± 0.02 


2.825 


1.64 ± 0.04 


1.64 ± 0.04 


4.64 ± 0.2 


1.89 ± 0.02 


2.875 


1.63 ± 0.04 


1.63 ± 0.04 


4.60 ± 0.2 


1.87 ± 0.02 


2.925 


1.63 ± 0.04 


1.62 ± 0.04 


4.56 ± 0.2 


1.85 ± 0.02 


2.975 


1.63 ± 0.04 


1.62 ± 0.04 


4.53 ± 0.2 


1.83 ± 0.01 


3.025 


1.63 ± 0.04 


1.61 ± 0.04 


4.49 ± 0.2 


1.81 ± 0.01 


3.075 


1.63 ± 0.04 


1.60 ± 0.04 


4.46 ± 0.2 


1.79 ± 0.01 


3.125 


1.62 ± 0.04 


1.59 ± 0.04 


4.43 ± 0.2 


1.78 ± 0.01 


3.175 


1.62 ± 0.04 


1.59 ± 0.03 


4.40 ± 0.2 


1.76 ± 0.01 


3.225 


1.62 ± 0.04 


1.58 ± 0.03 


4.36 ± 0.2 


1.74 ± 0.01 


3.275 


1.62 ± 0.04 


1.57 ± 0.03 


4.34 ± 0.2 


1.73 ± 0.01 


3.325 


1.62 ± 0.04 


1.57 ± 0.03 


4.31 ± 0.2 


1.71 ± 0.01 


3.375 


1.61 ± 0.04 


1.56 ± 0.03 


4.28 ± 0.2 


1.70 ± 0.01 


3.425 


1.61 ± 0.04 


1.55 ± 0.03 


4.25 ± 0.1 


1.68 ± 0.01 


3.475 


1.61 ± 0.04 


1.55 ± 0.03 


4.23 ± 0.1 


1.67 ± 0.01 


3.525 


1.61 ± 0.04 


1.54 ± 0.03 


4.20 ± 0.1 


1.66 ± 0.01 


3.575 


1.61 ± 0.04 


1.53 ± 0.03 


4.18 ± 0.1 


1.64 ± 0.01 


3.625 


1.61 ± 0.04 


1.53 ± 0.03 


4.15 ± 0.1 


1.63 ± 0.01 


3.675 


1.61 ± 0.04 


1.52 ± 0.03 


4.13 ± 0.1 


1.62 ± 0.01 


3.725 


1.60 ± 0.04 


1.52 ± 0.03 


4.11 ± 0.1 


1.61 ± 0.01 


3.775 


1.60 ± 0.04 


1.51 ± 0.03 


4.08 ± 0.1 


1.60 ± 0.01 


3.825 


1.60 ± 0.04 


1.51 ± 0.03 


4.06 ± 0.1 


1.59 ± 0.01 


3.875 


1.60 ± 0.04 


1.50 ± 0.03 


4.04 ± 0.1 


1.58 ± 0.01 


3.925 


1.60 ± 0.04 


1.50 ± 0.03 


4.02 ± 0.1 


1.56 ± 0.01 


3.975 


1.60 ± 0.04 


1.49 ± 0.03 


4.00 ± 0.1 


1.55 ± 0.01 
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Q [(breV/cj J 


r}2\}^ jxiU 


T}A\}^ jxiu 


jxiu 


rj8{(^ jxiu 




l.DU ± u.uo 


1 ACk 4- n OQ 

i.4y ± u.uo 


Q QS3 _!_ 1 

o.yo ± u.i 


1 f^A J- 01 
i.04 ± U.UI 


4.U / O 


1 KO -1- OQ 

i.oy ± u.uo 
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A 
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1 Al -\- n OQ 
1.4/ zt u.uo 


Q QO -1-0 1 
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i.OZ It U.Ui 




1 c^Q -1- n OQ 

I.oy ^ u.uo 


1 Al -\- OQ 
1.4/ ^ U.UO 


Q Q1 -1- 1 

o.yi ^ u.i 
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1.01 ^ U.UI 




1 CiQ -1- n OQ 
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1 AfK 4- n OQ 
1.4D ^ U.UO 
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1.4/ zt U.Ui 


/I 4 vr, 
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1 r,o -U n OQ 
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Q so -U 1 
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1 AO. 1 noi 

i.4D zt U.Ui 


i.r.)Z-.) 


1 ^si -1- n (V-* 
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g 1 -1- (1 1 
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4.0 1 
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±.0o di u.uo 
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1.44 ni u.uo 


Q 7Q -1- 1 

0. / y i u.i 
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1.40 ni U.UI 


4.DZ0 


l.Oo 3: u.uo 
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0. / ± U.i 
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Q -1-0 1 
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i.Oo It U.Ui 
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FIG. 12: Results ol the twist analysis. The squares represent the Nachtman moments obtained in this analysis. The solid line 
is the fit to the moments using Eg. 1201 with the parameters listed in Table ITTTl The twist-2, twist-4, twist-6 and higher twist 
(HT) contributions to the fit are indicated. The twist-2 contribution was calculated using Eq. 1211 



model dependent and we have checked these before ap- 
plying them to the data. Unfortunately, the Hall C data 
do not contain the quasi-elastic peak, therefore we had 
to compare to the previous SLAC and DESY measure- 
ments from Refs. |a^SIE3- Some of these data are not 
corrected for the radiative corrections, so we included 
radiative corrections in the model calculations. An ex- 
ample of the comparison of the quasi-elastic cross section 
model to the data, shown in Fig. 1141 indicates an overall 
few percent agreement and a particularly good match on 
the low W"^ side of the peak, which is important for the 
determination of the higher moments. 

At large values, the parameterization from Ref. 
fails to reproduce the data at high x. To correct the 
parameterization in this kinematic domain we switched 



to an older version of the fit reported in Ref. [53| for 
X > 0.75. However, in order to match the two fits we 
had to replace the x variable in the parameterization from 
Ref. [13 with x' =x- 0.9(x - 0.75)2. 

APPENDIX B: FIT OF THE RATIO R = olIot 

The ratio R{x, Q^) for the proton is well established 
in the DIS region and can be fairly well described by the 
SLAC fit from Ref. [s^. However, until recently, exper- 
imental data in the resonance region were missing. The 
data from Hall C published in Ref. [S^ cover the entire 
resonance region and extend down to very low values. 
It was shown by the HERMES Collaboration in Ref. 113 
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factor: 



q' (GeV/c)^ 




exp 



with = 0.8 (GeV/c)2, Cq = 0.729, Bq = 2.14, 
0.165 GeV, Bw = 0.383 and where Wth ^ M + is 
the value of the invariant mass at the pion threshold. 
Parameters listed above were obtained from the fit to 
experimental data on R taken from Refs. [sqIIgH. 

In this way the resonance structures, clearly seen on 
the proton, were averaged out to a mean curve. The dif- 
ference between these two models gave us an estimate of 
the systematic uncertainties of the ratio R, which turned 
out to be very small. 



FIG. 13: Ratio of the total higher twist (see Eq. |^ to the 
leading twist given in Ea. l21l with its systematic uncertainties. 
Stars - M2; squares - M4; triangles - Me; circles - Mg. 
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FIG. 14: Quasi-elastic cross section calculations in 
/ib/(GeV str) from Ref. [i^ compared to the experimental 
data from Ref. H at S = 9.766 GeV and 6 = 10°. The 
uncertainties are statistical only. 



that in DIS the ratio R does not depend on the nuclear 
mass number A. We take this assumption to be valid 
also in the resonance region. But since smearing effects 
of Fermi motion were expected to change both the F2 and 
Fl structure functions, we performed a "smooth" param- 
eterization of the measured ratio R. This smooth param- 
eterization is based on the fit from Ref. [s^, modified at 
low values and large x by means of a multiplicative 




0.9 1 
W in GeV 



FIG. 15: Deuteron ratio R — aL joT in the quasi-elastic region 
as a function of W at — 3.25 (GeV/c)^. The points are 
from Ref. [6^ and the curve represents the parameterization 
described in text. The minimum ai W = 1.07 GeV is due to 
the pion electroproduction threshold. 

The ratio R under the quasi-elastic peak is a separate 
issue. Because of the nature of the quasi-elastic peak, R 
is no longer independent of A and should therefore be 
treated within a nuclear model calculation. We used the 
model from Ref. 40^ , which treats separately the longitu- 
dinal and transverse nuclear response functions Rl and 
Rt to obtain the ratio R. The conventional ratio R can 
be calculated from those quantities as follows: 



R = 



2Q^ Rl 
Q2 + 



(B2) 



The deuteron quasi-elastic ratio R obtained from this 
model was compared to the data on the ratio R for the 
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deuteron ^ (see Fig. [TB|) and other nuclei (deuteron 
data on R in the quasi-elastic region are scarce) . Further- 
more, the calculations were compared to the sum of the 
proton and neutron form factors, which simply implies: 



R = 



f<2 



(B3) 



where r = Q'^/AM'^ and Gij, Gm are sums of the known 
Sachs form-factors of the proton and neutron. Since the 



number of protons and neutrons is different in different 
nuclei we did not expect to have the same ratio R for all 
of them. At the same time the x-shape of the ratio is 
very similar from nucleus to nucleus. In the low re- 
gion, where precise data exist, the calculations reproduce 
the x-shape of the data reasonably well. The systematic 
uncertainty was estimated as a difference between the 
model calculation and the result of the naive form-factor 
sum given by Eq. IB3I 
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